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Abstract 

We study the energy spectrum of the two-electron spherical parabolic quan- 
tum dot using the exact Schrodinger, the Hartree-Fock, and the Kohn-Sham 
equations. The results obtained by applying the shifted-l/N method are com- 
pared with those obtained by using an accurate numerical technique, show- 
ing that the relative error is reasonably small, although the first method 
consistently underestimates the correct values. The approximate ground- 
state Hartree-Fock and local-density Kohn-Sham energies, estimated using 
the shifted-l/N method, are compared with accurate numerical self-consistent 
solutions. We make some perturbative analyses of the exact energy in terms of 
the confinement strength, and we propose some interpolation formulae. Sim- 
ilar analysis is made for both mean-field approximations and interpolation 
formulae are also proposed for these exchange-only ground-state cases. 

I. INTRODUCTION 

Progress in nanotechnology has allowed the development of small devices like quantum 
dots. The confinement potential can be safely approximated by a harmonic one which 
has boosted the study of quantum dots with parabolic confinement during the last years 
(see e.g. Refs. 3 and 4, and references therein). The presence of many interacting electrons 
render the computation of the electronic states and properties a very complicated many- 
body problem. The first non-trivial exactly solvable problem of many-electrons is the one of 
two electrons confined in a parabolic potential, which made it a very attractive workbench 
for testing all kind of approximations (see e.g. PH10[). 
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The Hamiltonian describing the pair of interacting electrons in a parabohc quantum dot 
in the effective masss approximation can be written as 



where m* is the effective mass, u the confinement strength, and e the dielectric constant of 
the host material and is the Laplacian operator. 

Equation (|l|) can be separated into a centre-of-mass and a relative motion component. 
Furthermore, due to the radial symmetry of the components, only those parts of the cor- 
responding Schrodinger equations have to be solved. The centre-of-mass part will give the 
well-known harmonic oscillator problem. In three dimensions, the relative motion part may 
admit exact solutions for special choices of the parameters (see e.g. [|]). For the two di- 
mensional case, similar separation and substitution can be made, and again there are no 
general solutions expressible in terms of special functions. Nevertheless, in Ref. and 



more recently in Ref. it was shown that there exist analytic solutions for special choices 
of the confinement constant. 



Many-body effects are usually divided into exchange and correlation components [[14 
Exchange-only effects are considered in Hartree-Fock (HE) and differently, in Kohn-Sham 
(KS) approaches (although in KS, the correlation effects can be included), which typically 
amounts for around 10% of the total energy. Correlation is about one order of magnitude 
smaller. Nevertheless, it has been shown that in two dimensions for two electrons in a 



harmonic field, correlation may play a bigger role specially for singlet states (see e.g. ||15|| ). 

The shifted-l/N method [1^,1^ has been applied to various condensed matter problems. 
Also, the two-dimensional relative motion Schrodinger problem have been solved using this 
technique in Refs. P,p!8|. 

The article has been structured as follows: in Sect II we describe the Schrodinger, 
Hartree-Fock, and Kohn-Sham approaches we use. In Sect. Ill we solve the exact, HE, 
and KS-LDA equations using an accurate numerical technique and the shifted-l/N method. 
We also apply perturbation theory up to first order in both limits of confinement for the 
exact and mean field cases, and we propose some interpolation formulae for the energy. We 
discuss the accuracy of the mean field approaches, and of the shifted-l/N method for the 
present case. 



II. METHOD 

Throughout the paper the units of energy will be given in terms of the effective Rydberg 
constant TZ* = fi^ / {2m* a*"^) and the effective Bohr radius a* = ^^e/m*e^, respectively. 

In centre-of-mass and relative coordinates and measuring in reduced units the Hamilto- 
nian reads 




where we have chosen the centre-of-mass R = (n + r2)/2 and the relative coordinate f = 
(fi - f2). 
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The separation leads to the harmonic oscillator problem for the centre-of-mass coordi- 
nate, with energy 

E^L = ^{2N + L + 3/2) . (3) 

and eigenfunctions 

^7VLAf(i?) =Ar^Lexp(-7i?2)(27)V2^i-^^+i/2^27i?2)y^j,,(^^,0^) . (4) 

For the relative coordinate equation the wavefunction can be separated into radial and 
angular components 

vI/(r) = ^F,^(^,0), (5) 

where Yira{9, (p) are the spherical harmonics which are eigenfunctions of the angular momen- 
tum operator L^, and with eigenvalues m and /. This separation makes the corresponding 
radial Schrodinger equation the following second order ordinary differential equation 



1 o o 1 



dr^ 



+ /(/ + l)- + -7V + --E 



u{r) = . (6) 



It is well known that exact solutions of equation (|^) cannot be expressed in a closed form 
in terms of special functions. There are analytic expressions for the energy for particular 
values of 7 and / as it was pointed out in Refs. P,p!T|, among others. 

The electrons should satisfy the Fermi-Dirac statistics which means in this case that for 
singlet states [s = 0) the spatial part of the wavefunction should be antisymmetric and for 
triplet states (s = 1) symmetric. As the centre-of-mass coordinate remains the same after 
exchanging to electrons, the antisymmetry requirement will be in the relative part. Because 
of the separation in radial and angular components of the relative-coordinate wavefunction 
it will mean that singlet states are associated with odd / and triplet states with / even, 
respectively. 

It is interesting to compare the results of exact calculations with independent-electron 
models like Hartree-Fock (HF), and Kohn-Sham (KS) in order to assess the relative 
importance of many-body effects like exchange and correlation, and also to evaluate the 
performance of the local-density approximation. For two paired electrons the electronic 
density is phf = 2|(/)hfP, where (pHF is the orbital, and the exchange potential is equal to 
half of the Hartree one with opposite sign. The HF equation can be written as 

-^V^ + v{r) + ^vhIphf] 4>hf = £hf 4>hf , (7) 
f (r) = |7^r^ and Ehf is the HF orbital energy. The total HF energy is written as 

Ehf = '2b HF " ^ / drpHF vhIphf] , (8) 
where vh is the Hartree potential given by 
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dr 



(9) 



The KS equation can be written as 
1 



^V^ + t;(r) + vh[pks\ + v^{pKs) 



= £ks <Pks 



(10) 



Sks and (pxs are the KS orbital energy and eigenfunction, respectively, and again pxs = 
2|0K5p- We take here the local-density approximation for which Vx{p) = ^c^p^^^, with 



3 (3 



1/3 



see e.g. [0]). The total KS energy is given by 
Eks = '^£ks -^J drpKs vh[pks] ~ ^ / dfpxs VxipKs) ■ 



(11) 



Eqs. and (|T0D have asymptotae controlled by the harmonic potential, so the asymptotic 
density looks like 



Pa{r) oc exp(— 7r^ 



(12) 



III. RESULTS AND DISCUSSION 

First we make an analysis of the solution of the exact case using the 1 /N approximation. 
For details in the derivation of the formulae related to the shifted-l/N method, we refer to 
the literature (see e.g. [p|,|l6|-[T8[]). 



In order to establish the accuracy of the results obtained by the application of the shifted 
method we compare them with those obtained by using the Schwartz's numeric method 
|T^. The method is based in a numerical approximation of functions on a mesh and gives 
very accurate results pO| , pl| . There are only empirical estimates of the error which 
turns out to be exponentially decaying with the number of points given the mesh step. The 
interpolation function is chosen as 

/(r) = E/»^ , ^^'\ ^ (13) 

where 

u{r) = ?,m['K{r/hf'\ (14) 

Here is a zero of u{r), is a zero of its derivative, and h is the step of the mesh which 
turns out to be quadratically spaced. The choice of the step h was made after estimating the 
characteristic length of the effective potential, and then multiplying the obtained estimate 
by five and dividing it by the square of the number of points in the mesh, usually around 
300. This guarantees that h is minimal for a given uj, and also that the function value at 
the last mesh point is practically zero. 

In figure 1 we show the behavior of the error of the energy of few lowest eigenvalues 
calculated using the shifted-l/N method, compared with the accurate results obtained using 
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the Schwartz's method. We plotted the relative error defined as 5E = Eapprox/ E exact — 1, 
as function of the reduced variable 7' = 7/(7 + 1). The Os, Op, and Oc? states are the 
three lowest energy states of the relative motion with angular momentum / = 0, 1, and 2, 
respectively. It can be seen that the error remains bounded in 0.5% for the first eigenvalue, 
in 0.075% for the second one and 0.022% for the third one. It is noticeable that the method 
always underestimates the correct values of the energy, and that the error decreases with the 
increase of the angular quantum number I, as it should be expected. Furthermore, the error 
has some maximum between the two limiting cases, after which it decreases as expected, 
since the shifted-l/N method reproduces exactly the oscillator case. 

Figure 2 shows the behavior of the total energy (the sum of the centre-of-mass and 
relative motion energies) as a function of 7', as calculated using the shifted-l/N method. 
The first two symbols are the indices of the centre-of-mass component of the energy, and the 
last two correspond to the relative motion. The lowest six states are depicted. The inset is 
a magnification for 7 ' from 0.01 to 0.15. Here, an apparent linear behavior can be observed, 
but a more careful analysis indicates that, for instance, for the ground-state first six-eight 
points, the effective power in terms of 7 ' is about 0.8. In the last few points the curve can 
be fitted well with a function of the type 7 '/(I — 7') which is the expected behavior for 
electrons in a harmonic field. The relative error is not shown, but it is estimated in roughly 
one half of the one shown in figure 1, since the centre-of-mass energy can be calculated 
exactly and it is typically of the size of the relative motion energy. 

Although we are able to solve the Schrodinger equation very accurately for this system, 
perturbation analysis can give some more insight of the behavior of the electrons under weak 
and strong confinement. For weak confinement (7 — >■ 0) we have that the kinetic energy 
term can be neglected (see e.g. [13,22|), which corresponds to the strong interaction limit 



(Wigner crystallization). Then, the energy is approximately taken as the minimum of the 
effective potential (this is the zeroth order approximation to the energy), 

V^(r)=/(/ + l)i + ^7V + -, (15) 

the minimum is reached for Tq that satisfies the equation 

7^^^ - 2rl - + 1) = . (16) 

For Z = Tq = (2/7^)^/^ and for / ^ 7/2 Tq ^ {21 /'jY^'^. The minimum of the potential is 
then 

f^o = U^ + l)4i + ^M^ + 3^ (17) 

which is 2-1/33/2^2/3 for / = and 7(2/ + l)/2 + -^^/^/{2iy/^ for large enough values of /. 

The next order in the approximation is to get the effective frequency = 
irf2v5(r)/dr2| 



, 1 , 1 31(1 + 1) , , 



'0 



which is 7o = 3/47^ and 7^ = 7^ + 3/47^// + 2^/^/4(7//)^/^ for large /. Now we can estimate 
the energy levels of this effective harmonic field. So the total energy in the weak confinement 
limit is 
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Enl ^ /(/ + l)^! + 77^^' + \ + 2^i{n + 1/2) , (19) 



where n = 0, 1, ... For I = 



31/2 

For / large enough 



EnO ^ ^2-V3^2/3 + 3^2^^^ ^ ^^^^ ^ ^^^^ 



1 . , . /7\l/2 r 3^2 /^x3/2' 



E..-7(2/ + l)+ - +2 



1/2 



(n+1/2). (21) 



Equation (|20|) gives the explanation why the effective power of the ground-state energy for 
small 7' is approximately 0.8: it is between 2/3 and 1, the effective powers for the weak 
field limit for the relative motion and the centre-of-mass energies in terms of 7 ' ~ 7. 

In the strong confinement regime, the zeroth order approximation amounts to neglect 
the electron-electron interaction, so it corresponds to the oscillator's energy. Application 
of the first order of the perturbation theory for strong confinement (7 00) for the 
relative coordinate equation leads to 

Enl ^ E:r + {€T\r-'\KT) , (22) 

where E°f = 27(272 + / + 3/2) and substituting (p^f (which is similar to equation (|D ) into 
the above equation gives 

^1/2 

0'nl,2kO"n.l,2k' T{l + k + k' + 1)AC' (23) 

^ k,k'=0 

where ani,2k are the coefficients of the generalized Laguerre polynomials (of Eq.(|)) that 
satisfy the recursion 

k — n — 1 . 

^nl,2k = anl,2ik-l)^^^-^^Y/2) ^^^^ 

with Qnifi = 1- Furthermore, the normalization constant is given by 

n 

Ki^ = E ani,2kani,2k'Til + k + k' + 3/2) . (25) 

k,k'=0 

Based on the perturbative results of Eqs. (0), (|19|) and (|23| ) we propose the following 
interpolation formula 

,,-1 771(0) p(oo) 

7^ + 7 

where the superscripts (0) and {00) correspond to the zero and infinite confinement limits. 
For 7 — i> equation (p6D will return approximately the weak confinement limit and for 
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7 — ^ oo the strong confinement limit. For the zeroth order approximations equation 
leads to 

~ ^-1 + ^ ' ^"^'^ 

and for the first order pertubative results, it will be 

p(int), . l-'K + 27K^ + 1/2)] + 7[27(2n + / + 3/2) + 2-y^-f^/'A,,i] 

where Ani denotes the summation in the right hand side of equation (^). The interpolation 
scheme of equation (|27|) performs consistently bad, except for the very weak, and very strong 
fields. The error goes up to 41% for the ground-state, underestimating the correct values. 
The introduction of first order corrections, corresponding to the interpolation scheme of 
equation ( |2^ ) brings a dramatic improvement on the values: for the ground-state the relative 
error is never worse than 3.3% for the energy of the relative motion, which means that the 
relative error of the total energy is around 1.7%. 

In the case of the mean-field approximations like Hartree-Fock and Kohn-Sham, we have 
considered the paired-electron ground-state case, and since the confining potential is radially 
symmetric, the orbitals and the density are also radially symmetric. We implemented the 
shifted-l/N technique for the HF and KS equations. Here, a word about accuracy is needed: 
although the solution of the Schrodinger-like equations using Schwartz's method is very 
accurate, the estimate of the Hartree potential is not as accurate anymore, nevertheless, 
five to six figures are always guaranteed. The results are shown in table 1 in the first and 
third columns, indicated as HF — 1/N and KS — 1/N, respectively. Since the resulting 
wavefunction from the application of the 1/N has a complicated form which makes difficult 
a direct evaluation of the density or of the Hartree potential, we assumed that as initial 
guess the non-interacting density, which is also correct asymptotically. For comparison 
purposes, we have used again the Schwartz technique self-consistently to solve both the HF 
and KS equations. Numerical results are shown in the second and fourth columns of table 
1, indicated as HF — num and KS — num, respectively. Also for comparison purposes, we 
included the results of solving the full Schrodinger equation using the shifted-l/N method 
and the numeric solution, which are the last two columns indicated as Exact — 1/N and 
Exact — num, respectively. 

Table 1 shows very good agreement between the results from the accurate numerical 
result and the ones calculated with the 1/N method. We should not be too enthusiastic 
about the accuracy, since the remarkable agreement is probably a result of the compensation 
of errors from the calculation of the Hartree potential and energy (due to its simplicity) and 
the intrinsic error of the 1/N method, specially for the weak field case. It can also be seen 
that the relative accuracy improves from typically few times 10^^ to 10~^ — 10~^ with the 
increase of the strength of the field. Also the systematic difference between the HF and KS 
methods is reduced with stronger confinement. Both behaviors can be understood by taking 
into account that the confinement potential dominates over the decaying electron-electron 
interaction potential, and with the increase of the strength of confinement the problem 
becomes just a harmonic potential problem, for which the shifted-l/N is designed to give 
the exact energy, although the quality of the wavefunction is not too good. 
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In order to better understand the behavior of the pair of electrons in strong and weak 
confinement within these mean-field theories, we can recourse to pertubation analysis. In 
the strong confinement limit (7 — >■ 00) the system will behave basically as a pair of non- 
interacting electrons, and the electron-electron interaction (Hartree and exchange potentials) 
can be considered as a perturbation. The HF and KS orbitals become 

m = [I) exp(--7r2) , (29) 

so, the Hartree energy will be 2(27/7r)^/2 ^^^^ ^-^^ exchange 3c^2^^^{3-f/7rY^^ /A. Then 
the total HF energy will be 

4°°; = 37 + f ^V^' , (30) 



7r 



and the total KS energy 

Ei^^ = 37 + (2^/^ + h'/'S'/'c^^ (^) . (31) 

In the weak confinement limit (7 0), due to scaling arguments, we can neglect the 
kinetic energy, and we can assume constant density, at least within a certain radius R 
(this is asymptotically true for the HF approximation and arguably for the KS one). Now 
taking Eqs. ([Zp and (PD|), and using Poisson's equation, we get that phf ~ 7^/(27r) and 
pj^s 7^/ (47r). Then from normalization condition R = [3/(27rp)]^/^. The Hartree potential 
will take the form 

^H(r) = |-^, tfr<R, (32) 

and Vnir) = 2/r if r > R. The Hartree energy will be equal to 12/ (5i?). Substituting the 
above result into equation (H) we get that the Hartree-Fock energy will be 

E^^l ^ (37)^/^ . (33) 

Substituting the result for the Hartree energy into equation ([TI|) we find that in this limit 
the Kohn-Sham energy will be equal to 

E^l « -,=/=(6^V. + _|£_) . (34) 



Again, we can use interpolation schemes like the one of equation (|26|) , for the HF ground- 
state energy we will have 



In the Kohn-Sham case we can write down 

7^/3[63/^ + 2c./5/(47r)V3] + ^2^3^ ^ ^2^^ ^ 3 21/331/2,^4) (^/^)i/2] 
^Ks 17J ~ iTy ■ ^ ^ 
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Comparing the asymptotics of the exact equation for the ground-state in the strong 
confinement hmit with the asymptotics of the HF equation, we observe that they coincide. 
In the Kohn-Sham case the zeroth order is the same, but the coefficient of the first order 
perturbation is shghtly higher (1.62 compared to the exact a/2 1.41). In the low density 
limit (weak confinement), although the three asymptotes are proportional to 7^/^ in the 
zeroth order, the coefficients differ substantially, the exact is 1.19, the HF 2.08, and the KS 
3.17. This indicates a consistent overestimation of the energy by the independent-electron 
approximations, which is the expected behavior at least for the HF approximation. Our 
findings for this system for the weak and strong asymptotae, and the numerical results, 
suggest that KS-LDA energies are always higher than the HF ones, which is consistent 



with numerical experience on atoms and molecules flj]. From the results shown in table 



1 we can estimate the correlation energies for the ground-state for different confinement 
strengths: for 7' = 0.1, the correlation energy is about 5.6% of the total energy, meanwhile 
for 7' = 0.9 (strong confinement) it is about 0.17% (for the lightest many-electron atom. 
Helium, it is about 1.4%, and for Argon with Z = 18, it is only 0.14%). This adds evidence 
to the suggestion that for harmonic fields correlation effects are more important, also in 
three dimensions, especially for weak confinement. 



IV. CONCLUDING REMARKS 



In the present article we have calculated the energy spectrum of a two-electron spherical 
quantum dot for the few lowest states, using the shifted-l/N method and the very accurate 
numerical Schwartz's method. From the comparison of the numerical results we could assess 
the quality of the shifted-l/N method, which consistently underestimates the correct values, 
although the error is rather small, and it decreases with the increase of the relative angular 
momentum. We have also applied perturbation theory up to first order in both limits 
of confinement of the electron, and we have proposed some interpolation formulae for the 
energy. Inclusion of first order perturbation allowed to construct an interpolation expression 
that performs reasonably well. We also solved the mean-field Hartree-Fock and local-density 
Kohn-Sham problems for the ground-state. Using the shifted-l/N method we got from 
reasonable to high accuracy already in the first iteration, compared to the self-consistent 
numeric solution using the Schwartz's method. We made an analysis of the strong and weak 
confinement limits, and proposed interpolation formulae for both the Hartree-Fock and 
Kohn-Sham ground-state energies. It was shown that the correlation energy is relatively big 
for this systems, especially for weak confinement. 
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FIGURES 




FIG. 1. Relative error [SE = {Eijj^ / E^xact ~ 1) x 100) of the energy of the relative motion for 
the lowest energy states for Z = 0, 1, and 2 calculated with the 1/N method, as a function of the 
reduced frequency 7' = 7/(7 + l). 



11 



25 



20 



J.O 



15 



0.0 




0.0 



0.1 



Q ^ if- 

Hi 



f I c? 

9 i i 

/ J 



6 _l 



(OS,Os) + 

(OS,Op) X 

(OS,Od) * 

(OP,Os) □ 

(OP.Op) 1 

(OP.Od) 



10 



0.0 



0.2 



0.4 



0.6 



0.8 



1.0 



FIG. 2. Total energy as a function of 7 ' = 7/(7 + 1) for the few lowest energy states. 05, and 
OP depict the lowest energy states for the centre-of-mass motion for L = 0, and 1, respectively, 
and Os,Op, and Od the lowest energy states for the relative motion for I = 0,1, and 2, respectively. 
The inset shows a magnification for small 7'. 
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TABLES 
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w p ^ / AT 
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H/XdCZ — L 1 Is 


h/XCLct — nuTfi 


n 1 

U. i 


U.OOOD 


U.O / Oct 


u.oyou 


U.OUoZ 








l.llDo 






1 1 749 
1.1/ 


l.Uc50c5 


1 n»Q9fi 
i.uoyzo 


0.3 


1.7758 


1.7826 


1.8408 


1.8503 


1.7398 


1.74478 


0.4 


2.6200 


2.6255 


2.7029 


2.7118 


2.5791 


2.58569 


0.5 


3.7673 


3.7717 


3.8711 


3.8791 


3.7217 


3.73012 


0.6 


5.4477 


5.4508 


5.5775 


5.5842 


5.3972 


5.40775 


0.7 


8.1906 


8.1922 


8.3558 


8.3608 


8.1345 


8.14778 


0.8 


13.5693 


13.5693 


13.7902 


13.7928 


13.5057 


13.5232 


0.9 


29.3703 


29.3G79 


29.7094 


29.7082 


29.2930 


29.3194 



TABLE I. Three-dimensional two-electron quantum dot total energy, using the Hartree-Fock 
(HF), Kohn-Sham (KS), and exact Schrodinger (Exact) equations for selected values of the 
reduced confinement constant (7'), calculated using the shifted-l/N (1/A), and Schwartz's numeric 
(num) methods. 
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